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SUMMARY 

The neutral g-modes of a degenerate fluid at zero temperature are analy- 
zed. The g-modes of a degenerate fluid at finite but small temperatures are 
then expanded in terms of those of the zero temperature fluid. For non-relati- 
vistic degenerate fluids it is found that (a) the g-eigenvalues are proportional 
to T \i* y.\ where T is the internal temperature of the fluid, u and p. are 

e i “ \ 

the mean molecular weights of electrons and ions, respectively; (b) the ion 
pressure is solely responsible for driving the g-modes. For white dwarfs of 
about a solar mass,the periods of the g-osci 1 lations are in the range of a few 
hundred of seconds. 

I. INTRODUCTION 

It has been suggested that the short period oscillations observed in the 
cataclysmic variables are connected with the g-modes of the white dwarf component 
of these objects (e.g. Warner and Robinson, 1972 ; Osaki and Hansen, 1973 ). The 
thrust of the argument behind this suggestion is that the observed periodicities, 
of the order of several tens of seconds, are too long to be attributed to the 
p-modes. In connection with the longer periodicities of DA white dwarfs, of the 
order of several hundreds of seconds, some authors have also called upon the 
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g-modes. For example, Wolff (1977) has proposed that the Interaction of a slow 
rotation with the g-modes is capable of accounting for such periodicities. Most 
of these propositions have, at most, had partial acceptability by investigators 
in the field. Modifications to the theories have been proposed (e.g., Papaloizou 
and Pringle, 1978), and alternative explanations, not requiring the intervention 
of g-modes, have been put foreward (e.g., Bath 1973, and Bath et al., 197*0. 

In recent years, Baglin and Schatzman (1969) have calculated some g-fre- 
quenctes of white dwarfs with finite temperatures. Harper and Rose (1970). 

Osaki and Hansen (1973), and Brickhill (1975), each have analyzed the problem 
in varying details. Practically all previous works on non-radial oscillations 
of white dwarfs are attempts to integrate the differential equations governing 
the small displacements of the fluid. Some information on the eigenf requencies. 
and, to a lesser extent, on the eigenfunctions has been accumulated over the 
years. 

The authors wish to draw attention to an alternative approach to the 
problem: The operator generating small displacements of a fluid, including those 

of a degenerate structure, is self-adjoint (Chandrasekhar, 1964). The normal 
modes of the fluid belong to a Hilbert space. If there is access to a basis set 
for this space, then one can expand the actual eigendisplacements of the fluid in 
terms of this basis set, transform the differential equations of motions into a 
matrix equation, and obtain the expansion coefficients by variational calculations. 
This alternative route is followed in the present paper. It provides more infor- 
mation on the eigenvalues and eigenfunctions of the system, and gives deeper in- 
sight into their behavior and their dependence on the physical properties of the 
fluid. 
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(6) is an exact solution of Equations (1) corresponding to e = 0. This completes 
the demonstration of the existence of neutral modes in the zero temperature 
fluid, the neutral g-modes. 


IV. THE BASIS-SET 


The neutral state of the zero temperature fluid is infinitely degenerate. 

Any displacemertt vector ^ satisfying Equation (6), but otherwise arbitrary, is 
an eigenfunction of e s 0. Sobouti (paper l) has used this arbitrariness to 
propose a set of basis-vectors for the space of the g-modes of the neutral fluid. 
Thus, let t} t j = 1,2,...., have the following spherical harmonic expansion 

* ■ i - V H1 • i 

5 /i|r>(r) Y m 1 x J (r) - 3Y & 1 X J , (r) __L_ _L ) (11a) 

£ : ' r 2 ’ &U+1 ) r 30 ’ &U+1 ) r sin©. 3$ 

where x 1 = / dr • 

Substitution of Equation (lla) in Equation (6) gives 

x J ' = * j ' + f * j . <" b) 

One of the two functions V and X^ could be chosen arbitrarily. Equation (lib) 
can then be 'used to obtain the other one. One now appreciates the immense simpli 
fication that Equation (6), or equivalently Equation (lib), brings into the 
problem. These equations reduce the task of specifying a vector to the deter- 
mination of only one scalar function, ^ , say. The following expression for 


is used in the numerical calculations of this paper: 

j 3,-2 £+2j -2 5 _ i o 

pp p r » J - 1,2,.... 


(lie) 


We note the following: (a) On the surface of star t|» vanishes, while x', the 

non-radial component of the displacement, remains finite, (b) The exponent 
of r is chosen such that V.£ behaves as r £ near the origin and is an even or 
odd polynomial depending on whether l is even or odd; these properties are shown 
to be required by Hurley, Roberts and Wright (1966). (c) The polynomials 
in r are helpful in acheiving the completeness of the proposed vectors of 
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II. THE MODEL AND THE EQUATIONS OF MOTION 

A fluid consisting of non-relat ivi st ic degenerate electrons and non- 
degenerate ions will be considered. It will be assumed that the temperature 
is constant throughout the fluid, except in a thin envelope. In this envelope, 
the temperature will decrease to zero in a manner determined by whether the 
envelope is in convective or in radiative equilibrium. 

i) Equations of Motion: Adiabatic Lagrangian displacements of the fluid, 

1 /2 

?(r) exp(i e t) , are governed by the following equations: 


#</§ = ep£, ( la) 

where 

to/Z = V6p - 1 SpVp - pV6fi, ( lb ) 

Sp = -pV.E - Vp.E, (2a) 

Sp " ~ 'fp^'ad P M - 5>-£. (2b) 

2 

V = -ifiTGSp, (2cl 


and p, p and n are the density, pressure and gravitational potential of the 
fluid in hydrostatic equilibrium, respectively. 

ii) The Equation of State: The pressure of the fluid is the sum of the 

partial pressures of the degenerate free electronsand the non-degenerate ions. 
Thus , 

p(p.T) = P e (p. T ) + Pj(pJ)- (3a) 

No correction for charge separation is allowed in Equation (3a). The number 
density of the free electrons and the positive ions is assumed to be proport- 
ional to each other , so as to maintain macroscopic charge neutrality. There- 
fore, both p^ and p. are given in terms of the same density p. At temperatures 
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T, much below the Fermi temperature T p of electrons, temperature corrections 
In p e are of the order of (T/Tp) 2 . In contrast, the ions treated as a perfect 
gas have their pressure p ? proportional to the temperature. Thus, up to and 
including the first order terms in T/Tp, Equation (,3a) gives 


p(p,t) » p 6 (p) + ~hp t a nr pT F fi + if ’ Tp ' 


2/3 


(3b) 


where the first term on the right hand side is the electron pressure at zero 
temperature, y^ and y. are the mean molecular weights of electrons and ions, 
k and H are the Boltzmann constant and the mass of hydrogen atom, respective- 
ly. As a simplifying assumption, y g and y. will be considered to be constant 
throughout the star and also in the course of any thermodynamic process which 
may take place in it. 

iii) The Adiabatic Processes: The entropy of the fluid is the sum of the 

electronic and ionic contributions, s g and Sj, respectively. Thus, the change 
in the entropy is 


ds = ds„ + ds . . 

e l 

For the unit mass of ions, treated as a perfect gas, one has 


(4a) 


. >< (1*L _ dp_ x 

ds i ” y ; H 4 T p 


(4b) 


The entropy per unit mass of non-relativistic degenerate electrons is 

s = (iff 2 k/y e H) T/Tp (see for example, Landau and Lifshitz, 1959; Morse, 

I 965 ). From the last expression, and considering the fact that T p is propor- 
2/3 

tional to p , one obtains 


, „ 1 2 k T r dT _ 2 dp 1 

ds e 2 11 y e H Tp *- T 3 P > 


2 


IT 


3 


!il 

w e T F 


ds. . 
1 


(4c) 


One observes that (a) ds Is smaller than ds. by a factor of T/Tp, and (b) in 

the course of an adiabatic process (in which by definition ds * 0) both ds e and 
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ds. vanish simultaneously. These results simplify the forthcoming calculations 
considerably. Thus, one obtains 


(9P) 

Dp'ad 


- (SPb.) 

3P ! 


mi 

3P si 


„ 5 
3 


Pe 

P 


y£I 

p 


where y = 5/3 is ratio of the specific heats of the ion gas. 
III. THE NEUTRAL MODES OF THE ZERO TEMPERATURE FLUID 


(5) 


At zero temperature the g-modes of a degenerate fluid are neutral, that 
is, the g-eigenval ues become zero (e.g., Harper and Rose, 1970). We intend to 
analyze these neutral modes in some detail and subsequently use them as a basis 
set to expand the modes of the finite temperature model. The formalism is the 
same as that of Sobouti (1977a, henceforth referred to as paper l) for the g- 
and p-modes of ordinary fluids. 

Let us investigate if there is a displacement ?(r) ? 0 which is a solution 
of Equation (l) and for which the corresponding 5p is zero. This last condition, 
by Equation (2a), requires that 

6p = -PV.C - Vp.£ =0. (6) 

Substitution of Equation (6) in Equation (2b) gives 

6p = "(0ad " ( |?sJpZ*5» (7) 


where we have utilized the fact that in hydrostatic equilibrium p is a func- 
tion of p alone and VP = (3p/3p) Vp. The expression (3p/3p) „ is the deriva- 
tive of p with respect to p as prevailing in the equilibrium structure. The 
subscript "st" is inserted to distinguish this derivative from the adiabatic 
derivative (Sp/Sp)^. From Equation (3b) one has 


(lR) 

v 3p'st 


& 
V 3p ' 


st 


y jH 


k ,£Tx 
y jH pv 3p ; st ’ 
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(8a) 



The first derivative on the right hand side of Equation (8a) will be calculated 
for the zero-temperature structure and is equal to j P e /P* w,th regards to 

(3T/Dp) st , let us define y' as follows: 

Ost ■ <*' - »?• <8b) 

We observe that in the isothermal core, y' = 1. In the envelope, where the 
temperature decreases outwards, one has the following information: (a) In any 

region in convective equilibrium, y'>y . (b) In any convectively neutral 

region, y' = y. (c) In any region in radiative equilibrium, y'<y. From the 
theory of g-modes in ordinary fluids one knows that in case (a), some unstable 
g-motions (or g-modes, if they develop into standing patterns) arise. These 
are, however, confined to the convective layers. In case (c) stable g-modes 
develop and again are confined to the radiative region. We shall assume that 
the envelope is convectively neutral, case (b). By this assumption the envelope 
will not contribute to the g-modes of the system and the role of the degenerate 
core will be singled out. In view of these considerations, Equation (8a) becomes 

(|£) = — + EL in the core, (8c) 

v 9p'st 3 P P 

= ESL + y' EL , In the envelope. (8d) 

3 P P 

Substitution of Equations (5) » (8c) and (8d) in Equation (7) gives 

Sp = -(y - 1)p. V.£ = -(y - 1) ~pr T PV.£, in the core, (9a) 

= 0, in the envelope (y = y0. (98) 

This completes the reduction of 6p. From Equations (2c) and (6) one has 

= 0 . 0 °) 

We now observe that at zero temperature, 5p,6p, and generated by the 
displacements prescribed by Equation (6), vanish identically. Thus £ of Equation 
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Equations ( 1 1 ) in the g-subspace of the normal modes of the neutral-fluid. The 
completeness of the proposed vectors, though not yet established theoretically, 
has, however, been born out by numerical computations of Sobouti (1977b) and of 
this paper. Variational calculations with the basis vectors of Equations (ll) 
enable one to isolate the g-modes of a fluid systematically and with satis- 
factory accuracies. 

We have only discussed the g-modes of the neutral fluid. A second basis 
set for the p-modes can be generated from the requirement that the p-set is 
orthogonal to the g-set. These two basis sets may then be used to expand the 
eigend isplacement vectors of any other fluid at finite temperatures (see paper 
I for the p-basis set and for further details). Silverman and Sobouti ( 1978 ), 
and Sobouti and Silverman ( 1978 ) have carried out such an expansion for ordinary 
fluids. Their analysis shows that in the limit of small departures from the 
neutral state (the limit of small T in the present problem), the g-modes, as 
given in terms of the vectors of Equations (ll), are independent of the p-basis 
vectors. This property stems from the fact that the g-states at T = 0 are 
degenerate. In the language of linear vector spaces, to a degenerate state 
there corresponds a subspace of the normal modes of the system. The effect of 
a small perturbation is to specify the principal directions of this subspace 
but it leaves the subspace unaltered. Thus, up to the order T/T p , the g-eigendis- 
placements can be expanded in terms of the basis set of Equation (ll) alone. 

No intervention of the p-basis set will be necessary. 

V. THE EQUATION GOVERNING THE g-MODES: 

£ » J = 1,2,.... be a sequence of the g-eigendisplacements of the 
finite temperature fluid. The Eulerian variations < 5 p , 6p and Stt, generated by 
these vectors are given by Equations (6), (9) and (10), respectively. 

Substitution of these in Equations (1) gives 
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tVJ*} = V5 j p = e J ‘p£ j 


The density p in Equation (12), and in Equation (9) for 6^p, is the density of 
the zero temperature fluid. This does not imply, however, that the density and 
the pressure are approximated by their values at T a 0. We have been able to 
show that Equation (12) is correct up to the terms of first order in T/Tp, The 
terms arising from 6p, 6p and 6J1 in Equations (1), on account of the fact that 
the density distribution at T 5 * 0 is different from that at T = 0, cancel out 
each other. 

Let have the following expansion in terms of the basis set of Equations 

( 11 ): 

£ j = ZC k Z kj ; k, j = 1,2 03) 

The expansion coefficients Z k- i can be considered as the elements of a matrix 


Z. Thus, 


Z = [Z kj ] ; k, j = 1,2,.... 


0*»a) 


Each column of Z is an eigenvector of the system. Let E be the diagonal matrix 
of the eigenvalues , 


E = [ e jj diagonal, j . , >2 

In connection with Equation (13), let us introduce the following matrices: 


(l*b) 


W kj = J C k *.7(/C J *dv, (15a) 

S kj = Jp^. £ j dv. (15b) 

We are now ready to convert the differential Equation (12) into an equivalent 
matrix equation. Let us substitute Equation (13) in Equation (12), premultiply 
the resulting equation by C (say), and integrate over the volume of the 
fluid. We obtain the (i,j)th element of the following matrix Equation: 


WZ 


SZE. 
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The matrix Z simultaneously diagonalizes W and 5. A normalization condition 
on the eigendi splacements C J , or on the eigenvectors Z, could be imposed such 
that S diagonalizes to a unit matrix I (see Silverman and Sobouti, 1978). Thus 


Z + SZ = I . 


(16b) 


Given the matrices W and S, the eigenvalue matrix E and the eigenvector matrix 
Z can be solved from Equations (16). 


The W-and S-Matrices: From Equations ( 1 5a) , ( 1 2) , and (9), aft 
tegration by parts, one obtains 


er an in- 


=j(y - y') 


p.V. ? " V.£ J dv, 


(17a) 


where y 1 is defined by Equation (8b). The value of y* is equal to unity in 
an isothermal core; and is equal to y, the ratio of specific heats of ions, in 
a convectively neutral envelope. Substitution of Equation (11) in Equation 
(17a) and an integration over the solid angle' gives 



(17b) 


where R is the radius of the star. The assumption of a convectively neutral 
envelope confines the domain of integration to the core of the star, (i) The 
thickness of the envelope, however, is proportional to T and to R/M. 

The latter is, in turn, proportional to T" 1 . Therefore, the ratio of the 
thickness of the envelope to the radius of the star is of the order 
of T/Tp c . (ii) The density, temperature, and pressure tend to zero 

at the surface. These two factors make. the contribution of the envelope to the 
W-matrix of the second order in T/T^. Therefore, extending the domain of 
integration over the whole star does not introduce first order errors. This 
completes the reduction of the W-matrix. For the S-matrix, after an integration 
over the angles, one gets 
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1 k* j S . 

Tti+n x x J dr - 


08 ) 




The symmetry and positive definiteness of the W-and S-matrices are manifest from 
Equations (17) and ( 1 8 ) . 


VI. SOME PROPERTIES OF THE EIGENVALUES AND EIGENVECTORS 


i) The Sign and the Asymptotic Behavior of the Eigenvalues: From Equations 

(12) and (9) one can readily write down the following integral expression for : 


e J = (y - 1) •— \ pV.£ J “v.£ J "dv / ( p^'.^dv. 09) 

Either from Equation ( 1 9) or from Equations (17) and (18) for the W-and S-matrices, 
one can draw the following conclusions: 

a) The numerator and the denominator in Equation (19) are both positive 
for any , or, the W-and S-matrices are symmetric and positive definite. There- 

r * 

fore, all eigenvalues are positive. That is, the g-modes of a degenerate 
fluid are all stable. This, of course, is not surprising. The g-modes are 
prototypes of convective motions. An isothermal core is highly subadiabatic and 
is stable against convective motions. 

b) In Equation (19), as the mode-number j increases, the numerator becomes 
progressively smaller than the denominator. This can be seen from Equations 
(17b) and ( 18 ), where the integrations over the solid angles have been carried 
out. This asymptotic behavior is well estabilished for the g-modes of ordinary 
fluids (e.g. ,• Ledoux and Walraven 1958). 

ii) The Unit and the Order of Magnitude of the Eigenvalues: In Equation 

( 19 ), let us use a dimensionless radius variable x = r/R. Thus, one obtains 


= J - w- 'VMr- 


1 


R Fc 


(20a) 
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6 1 m j p !x\£ J V S Jd3 * / 1 e£J*- S J d3x, 


(20b) 


where Tp c is the Fermi temperature at the center. |n Equation (20a) the factors 
(T/T Fc ) and gJ are dimensionless. The first of these factors depends on the 

model under study. The second factor is calculated for a zero temperature star 
and is given in Tab1esl,2, and 3. The remaining factor in Equation (20a) has 

physical dimensions. On expressing T Fc and R in terms of the central density 


p c (see Chandrasekhar, 1 939-) one gets 

V, 


, AirGp 

cJ - 

n 

l 


( Y - ,) T 

P i T Fc 


j 


(21a) 


where is the Emden radius of the polytrope 1.5 (we note that the structure of a 
non-relativistic degenerate fluid at zero temperature is that of the polytrope 

1.5). The p-eigenvalues of the fluid are of the order of i»TrGp / n 2 . There- 

fore, the g-eigenvalues of the white dwarfs are smaller than the p-eigenvalues 

by a factor of T/T pc . Upon expressing p c in terms of the total mass of the star, 
one obtains an alternative expression for e-*" (see Chandrasekhar, 1939), 


(21b) 


J = 128tt 2 G 4 -6 i d0 x-2 / , * 5 u M 9 _ 1 

e 12^ q *1 1 \ j vT 'i vY " Du \ 2 T c. J 

125 K 3 1 dn 1 * pr'He ) r ’ 

Fc 

where (d0/d n ) j is the surface value of the derivative of the Emden temperature, 
and k is defined by the relation p = K(p/u e )5/3. The oscillation period 
P J = 2ir//eJ, in the cgs units, is given by 


PJ (sec) = 69.0(y - l) -1 / 2 y-5/2 (^Mjl/2 Mo j 1/2 , j.-l/2 

e Pe "FT ( j • ' ( ' 

Por y - 5/3 » y e = 2, and p. = 4 , one has 


(22a) 


PJ tsec) = 21.1 te J)-'/2 (22b) 

iii) Behavior of the Eigendisplacements at the Center and at the 
Surface: Once the eigenvector matrix Z is calculated from Equation (.16), the 
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eigendisplacement vectors £ J ' , can be obtained by Equations (13) and (ll). We 

note the following properties of £•* : 

a) At the center V,£ tends to zero as r^. See the remarks following 

yv ^ 

Equation (11c). 

b) At the surface, the radial component of £ tends to zero as r-R. The 

non-radial component and remain finite. These requirements from the g- 

modes of any fluid are analyzed in paper I. 

VII. NUMERICAL RESULTS 

Equation (l 6) is solved in a Rayleigh-Ri tz approximation, that is, by 
approximating the infinite matrices of Equation (13)-(l6) b Y finite matrices. 

The matrix size was varied from one by one to seven by seven. Values of 
l = 1,2,3,i»,5 and 6 were considered. For £ = 1 and 2 the dimensionless eigenvalues 
Equation (21b), and the eigenvectors Z at various Rayleigh-Ritz approximations 
are given in Tables 1 and 2. The eigenvalues are displayed in lines marked by 
an asterisk. The column below each eigenvalue is the corresponding eigenvector. 

The eigenvectors are normalized according to Equation (16b). For £ = 3,^,5, and 
6 the dimensionless eigenvalues are given in Table 3. The convergence of eigen- 
values is satisfactory. The same ansatz, Equation (11c), for non-degenerate 
fluids, however, gives much faster convergence (Sobouti, 1977b). Motivated by 
this observation a search for a more suitable ansatz for degenerate fluicTs is 
being made. 

VIII. CONCLUDING REMARKS: 

The salient features of our model can now be summarized. The white dwarf 

core is taken to be an isothermal fluid, composed of ions and non-relat i vi st i c 

degenerate electron gas. In this formalism, corrections up to first order in 

T/T are included, where T is the temperature of the core and T is the Fermi 
Fc 
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Table 1. The eigenvalues and the eigenvectors of the g-modes of white 
dwarfs, corresponding to i. = 1. The eigenvalues (Equation 20b) are displayed 
in rows marked by asterisks. The column below each eigenvalue is the corres- 
ponding eigenvector. Computations are in the Rayleigh-Ri tz approximation 
using from one to seven variational parameters. The last two digits in the 
entries represent the exponents of 10. 


* 0.6619-01 

- 0 . 1690-01 
0. 4082+00 
-0.2958+01 
0.8779+01 
- 0 . 11 2^+02 
0.5136+01 


0.5163-01 0.7871-01 


-0.1879-01 0.7072-02 

0.5878+00 -0.2476-01 
-0.5740+01 -0. 1264+01 
0.2444+02 0. 1042+02 
-0.5036+02 -0.2968+02 
0.4927+02 0.3544+02 
- 0 . 1830+02 - 0 . 1508+02 


9 7 


* 0.5l8i»+00 

0.6600-02 


* 0.2654+00 0.6724+00 

-0.8996-02 0.1352-02 
0.3167-01 0.2471-01 

* 0.1684+00 0.3497+00 O.7312+OO 

0.1109-01 -0.2199-02 0.3254-02 
-0.8377-01 -0.3596-01 -0.3898-02 
O.II83+OO 0.1072+00 0.6105-01 


0.8629-01 

- 0 . 1501-01 

0.2673+00 

-0.1340+01 

0.2449+01 

-0.1452+01 

0. 1043+00 

0.5793-02 

0.1235-01 

-0.9435+00 

0.4760+01 

-0.8043+01 

0.4376+01 

0.1190+00 

-0.9869-02 

0.1642+00 

-0.1440+01 

0.7994+01 

-0.2277+02 

0.2926+02 

-0.1347+02 

9 5 


0.1172+00 

-0.1307-01 

0.1607+00 

-0.4936+00 

0.4151+00 

0.1461+00 

-0.4496-02 

-0.3400-01 

0.6051+00 

-O.1707+OI 

0 . 1283+01 

0.1673+00 

-0.8360-02 

0.1018+00 

-0.7592+00 

0.3378+01 

-0.6298+01 

0.3830+01 

0. 1842+00 

-0.3549-02 
- 0 . 6866 - 01 . 
0.1195+01 
- 0 . 6686+01 
0.1854+02 
-0.2444+02 
0 . 1185+02 


0.2166+00 

0.3303-02 

,0.4027-01 

-0.3102+00 

0.3772+00 

0.2500+00 

0 . 6698-02 
-0.5606-01 
0.3828+00 
- 0 . 12 Z 7+01 
0.1104+01 

0.2743+00 

-0.2535-02 

-0.5969-01 

0.6453+00 

-0.2605+01 

0.4901+01 

-O.3225+OI 

0.2924+00 

0.5469-02 
-0.6653-01 
0.8952+00 
-0,5338+01 
0. 1475+02 
-0.1953+02 
0 . 9812+01 

9 3 


0.3999+00 

-0.4963-02 

0.2508-01 

-0.1706+00 

0.2980+00 

0.4288+00 

-0. 1810-02 
-0.4056-01 
0.2621+00 
-0.7487+00 
0.7824+00 

0.4471+00 

-0.3873-02 

0.3168-01 

-0.4014+00 

0.1705+01 

- 0 . 3108+01 

0 . 2161+01 

0.4594+00 

-0.1605-02 

-0.4936-01 

0.5724+00 

-0.3339+01 

0.9375+01 

-0. 1237+02 
0 . 6348+01 

g 2 


0.7494+00 

0.1705-02 

0.1985-01 

-0.4368-01 

0.1223+00 

0.7554+00 

0.2531-02 

-0.1341-02 

0.9413-01 

-0.2043+00 

0.2482+00 

0.7580+00 

0.1855-02 

0.1848-01 

-0.8743-01 

0.4633+00 

-0.7988+00 

0.5804+00 

0.7593+00 

0.2360-02 

-0.1686-02 

0.1523+00 

-0.7685+00 

0.2229+01 

-0.2916+01 

0.1519+01 
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Table 2. The eigenvalues and the eigenvectors of the g-modes of white 
dwarfs, corresponding to 1-2. See the legend for Table. 1 for further details 


* 0.1199+01 

0.21A8-01 

* 0.6630+00 0.1567+01 

-o. 351*8-01 -0.1502-02 
0.9286-01 0.7693-01 

* 0.1*339+00 0.8694+00 0.1733+01 

0.5113-01 -0.921*1-03 0.1076-01 
-0.2908+00 -0.1505+00 -0.1*392-01 
0.3389+00 0.3096+00 0.1990+00 

* 0.3077+00 O.5522+OO 0.1002+01 0.1800+01 

0.6907-01 0.5869-02 -0.2000-01 -0.2757-03 

-0.61*27+00 0.2205+00 0.1333+00 0. 881 1*-01 
0.1635+01 -0.1083+01 -0.6853+00 -0.2668+00 
-0.1200+01 0.1080+01 0.8935+00 0.1»55l»+00 


* 0.2299+00 0.3798+00 O.6367+OO 0.1088+01 0.1830+01 

0.891*6-01 -0.1280-01 0.3036-01 0.l!*61»-02 0.8011*-02 
-0.1208+01 -0.2673+00 -0.273 z »+00 -0.2l8l»+00 -0.6347-01 
0.5007+01 0.21*99+01 0.1672+01 0.1198+01 0.531*5+00 

-0.7972+01 -0.5668+01 -0.1*258+01 -0.2937+01 -0.1 1 67+01 

0.1*257+01 0.3722+01 0.3237+01 0.21*99+01 0.1085+01 

* 0.1783+00 0.2765+00 0 . 1*334+00 0.6996+00 0.111*6+01 0.181*6+01 

0.11 2l*+00 0.2150-01 -0.1*239-01 0.8077-03 -O.I678-OI 0.2506-03 

-0.2053+01 0.2660+00 0.1*999+00 0.3583+00 0.221*8+00 0.1108+00 

0.1227+02 -0.1*656+01 -0.351*9+01 -0.2946+01 -0.2093+01 -0.7626+00 
-0.3160+02 0.1816+02 0.1 32l*+02 0.1036+02 0.7l»59+01 0.2920+01 

0.3630+02 -0.261»l»+02 -0.2105+02 -0.1683+02 -0.1200+02 -0.1*595+01 
-0.1523+02 0.1296+02 0.1136+02 0.9706+01 0.7239+01 0.2855+01 

* 0.11*08+00 0.2105+00 O.3129+OO 0.1*770+00 0.7500+00 0.1191+01 0.1857+01 

0.1377+00 0.3290-01 0.5618-01 -0.5112-02 -0.262l*-01 0.3197-02 0.7393-02 

-0.3255+01 0.1637+00 -0.851*8+00 -0.5043+00 0.1*21*6+00 -0.3363+00 -0.9939-01 

0.2617+02 -0.7301+01 0.7055+01 0.5960+01 -0.1*739+01 0.3396+01 0.1279+01 

-0.9649+02 0.1*41*3+02 -0.3412+02 -0.2857+02 0.2349+02 -0.1684+02 -0.6090+01 

0.1780+03 -0.1078+03 0.8399+02 0.6946+02 -0.5700+02 0.4117+02 0.1507+02 

-0.1594+03 0.1152+03 -0.9583+02 -0.8137+02 0.6736+02 -0.4870+02-0.1778+02 

0.5511+02 -0.4508+02 0.4019+02 0.3571+02 -0.3048+02 0.2246+02 0.8291+01 

g ? g 6 g 5 9* 9 3 ^2 9 i 
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Table 
ponding to 

z = 3 

0 . 2612+00 
Z = 4 


0.3997+00 
z - 5 

0. 5445+00 
£ = 6 

0 . 2131+00 

S 7 


3- The g-eigenvalues of white 
£=3,4,5, and 6. See the legend 


0.5477+00 
0.4135+00 0.6708+00 
0.3237+00 0.4941+00 0.7637+00 
0.3788+00 0.5557+00 O.8381+OO 


0.8231+00 

0.6266+00 0.1003+01 
0.4947+00 0.7447+00 0.1139+01 
0.5785+00 0.8341+00 0.1250+01 


0.1126+01 
0.8620+00 O.1366+OI 
0.6841+00 0.1021+01 0.1551+01 

0.7978+00 0.1151+01 0.1712+01 


0.1450+01 
0.1115+01 0.1755+01 
0.8802+00 0.1320+01 0.1998+01 
0.1089+01 0.1512+01 0.2174+01 

g 6 9 5 9 <. 


dwarfs (Equation 2 1 6) , corres- 
for Table 1 for further details. 


0.2005+01 
0.1146+01 0.2645+01 

0.7623+00 0.1496+01 0.2967+01 

0.9626+00 0.1733+01 0.3122+01 

0.1109+01 0. 1 896+01 0.3204+01 

0 . :1 221+0 1 0.2014+01 0. 3254+01 

0.1310+01 0.2105+01 0.3289+01 

0.2942+01 
0.1694+01 0.3916+01 

0.1136+01 0.2206+01 0.4444+01 
0.1427+01 0.2565+01 0.4724+01 
0.1643+01 0.2823+01 0.4888+01 

0.181 2+01 0.3017+01 0.4995+01 

0.1952+01 0.3173+01 0.5074+01 

0.4023+01 
0.2299+01 0.5395+01 

0.1546+01 0.2993+01 0.6174+01 

0.1934+01 0.3492+01 0.6614+01 

0.2229+01 0.3862+01 0.6887+01 
0.2463+01 0.4149+01 0.7074+01 

0.2669+01 0.4393+01 0.7221+01 

0.5260+01 
0.2960+01 0.7096+01 
0.1986+01 O.3856+OI 0.8166+01 
0.2481+01 0.4514+01 0.8798+01 
0.2861+01 0.5014+01 0.9205+01 
0.3176+01 0.5417+01 0.9500+01 

0.3359+01 0.5615+01 0.9625+01 

9 3 9 2 9 1 
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temperature at the center. 

From Equations (.21), it is seen that the eigenvalues for g-mode oscil- 
lations are proportional to the internal temperature T. We also note that the 
factor (y - 1) appears in Equation (,21a) because of the assumption of an 
isothermal core. If there is any decrease in the temperature outwards, this 
factor should be replaced by (y - y'), where y 1 ' s defined by Equation (8b). 
This replacement will result in a decrease in the eigenvalues. 

Again, from Equation (21b) we see that the eigenvalues are very sensitive 

.6 

functions of the mean moleculaij weight of the electrons, varying as y g . 
Dependence on y., the mean molecular weight of the ions, is not so pronounced. 
These features should be comparecl\with B r i ckh ill's remark that the periods of 
gravity oscillations of white dwarfs do not depend critically on the composition 
of the stars. 

The periods of g-modes are given by Equation (23), where the values of 
are of the order of unity (see Tables 1,2, and 3). Thus, the periods of 
oscillations come out to be of the order of a few hundred seconds. As the 
harmonic number £ increases, the periods decrease. These conclusions, as well 
as the linear dependence of our eigenvalues on the temperature are in agreement 
with similar results obtained, through an entirely different approach, by 
Papaloizou and Pringle. 

Let us also emphasise an important difference between the g-modes of a 
degenerate structure and those of an ordinary fluid. According to Cowling 
(19^1) > the pressure fluctuations associated with the g-modes of an ordinary 
fluid are less prominent than the corresponding density fluctuations. For 
degenerate fluids, the opposite appears to be the case. Equations (12) and 
(9a) show that the Eulerian changes in ionic pressure are mainly responsible 
for the g-modes. The effects of Eulerian changes in the density and in the 
electronic pressure remain unimportant (in the limit of non-relativistic 
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degeneracies). That the g-modes are primarily driven by the buoyancy forces 
remains valid for degenerate as well as for non-degenerate fluids. 

These calculations have been carried out for non-relativistic fluids. 

There are strong indications that in a relativistic case, the oscillation periods 
will decrease. 
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